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Abstract 

Recent photoabsorption measurements have revealed a rich fine structure 
in the cohective charge-density excitation spectrum of few-electron quantum 
dots in the presence of magnetic fields. We have performed systematic com- 
putational studies of the far-infrared density response of quantum dots, using 
time-dependent density-functional theory in the linear regime and treating 
the dots as two-dimensional disks. It turns out that the main characteris- 
tics observed in the experiment can be understood in terms of the electronic 
shell structure of the quantum dots. However, new features arise if a break- 
ing of the circular symmetry of the dots is allowed, leading to an improved 
description of the experimental results. 
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I. INTRODUCTION 



The electronic structure and dynamics of quantum dots have been a subject of intense 
study in recent time.l'i Frequently one makes the simplifying assumption that the electrons 
in quantum dots form a two-dimensional system confined by a parabolic potential of circular 
symmetry. The electronic ground state can then be found either through diagonalization 
schemes or using mean-field approaches such as Hartree theory or density-functional theory. 
The former treatment is computationally much more demanding and has so far only been 
used to describe few-electron systems, whereas the latter approach has been applied to 
electron numbers of order 100. In either case, however, one finds that the main features 
observed in the experiment (such as effects related to the electronic shell structure) are well 
reproduced. 

The model of a parabolic circular quantum dot is of course an idealization. It has already 
been recognized in the early 90 's that explaining certain experimental features requires 
including some deviation from parabolicity of the confining potential. A first study of this 
kind has been performed by Broido et al.p who considered the ground state and far-infrared 
(FIR) response of up to 30 electrons confined in a circular quantum disk of radius 100 nm. 
Their approach, as well as a similar one by Gudmundsson and Gerhards,@ was based on a 
Hartree plus RPA description of the electrons. Both studies showed that a correction to 
the confining potential towards the edge of the dot that makes it steeper than the parabolic 
potential leads to a blue shift of the collective dipole modes that increases with the number 
N of electrons in the dot. In addition, it was found in Ref. ^ that one has to include 
a deviation from circular symmetry in order to explain the anticrossing behavior in the 
FIR spectra observed in the experiments by Demel et a/.l These findings were supported 
by an exact diagonalization study of quantum dot helium by Pfannkuche and Gerhards.§ 
We finally mention a very recent study by Hirose and Wingreeni that uses spin-density- 
functional theory to describe the electronic shell structure and calculate addition energies 
in elliptic dots. 

The present work has been motivated by a recent series of experiments on self-assembled 
InAs quantum dots performed in the group of J. P. Kotthaus in Munich. The experiments 
use a combination of in situ capacitance spectroscopy and FIR absorption spectroscopy to 
probe the ground state and collective charge-density excitations as a function of the electron 
number per dot, 1 < < 6. For A^ = 1, 2 one finds the well-known two- mode spectrum 
described by the simple formula 



where c^o is the characteristic frequency associated with a parabolic confinement potential, 
and Uc = eB/m*c is the cyclotron frequency for particles of charge e and mass m* in a 
magnetic field B. However, for N > 2 the authors of Refs. detect a much richer FIR 
spectrum: the a;+ mode is observed to split up into three subpeaks (see Fig. 3 of Ref. |^). 

The authors argue within a simple picture of non-interacting particles that the observed 
splitting of the mode is caused by the non-parabolicity of the confining potential of 
the quantum dots. Since for A^ = 3,4,5 the p-shell is only partly filled, the system can 
perform transitions of the s ^ p and p ^ d type, which have slightly different energies, 



2 



in contrast with the strictly parabohc case. This simple explanation, however, can only 
account for a two-fold splitting of the 0;+ mode. It is then further argued that the observed 
third (somewhat weaker) signal is caused by effects related to electron-electron interaction. 

It is the purpose of the present paper to provide a theoretical explanation of the three- 
fold splitting of the ujj^ mode observed in Ref. Our approach is based on density-functional 
theory for the ground state and linear response in the presence of a magnetic field. It will 
turn out in the following that including electron-electron interaction effects alone is not 
sufficient. Instead, the presence of the third subpeak is explained by a combination of the 
non-parabolicity plus breaking of the circular symmetry of the confining potential. 

The paper is organized as follows: in Section 2, we shall introduce our model for the 
quantum dots, a two-dimensional disk, and we shall present the theoretical methods for 
describing the electronic ground state and dynamics. In Section 3, we then discuss our 
numerical results, and Section 4 gives our conclusions. 



II. FORMAL FRAMEWORK 
A. Ground-state calculation 

The electronic ground state of a system in the presence of magnetic fields is described 
by current-density-functional theory. The general formalism has been developed by Vignale 
and Rasolt,0 and to date there exist several applications to quantum dots.S'El Here, we 
want to generalize the formalism, which originally had been developed for circular quantum 
dots, to describe non-circular systems. 

In the following, we shall make our lives a little easier, especially in view of the linear 
response calculations later on, and neglect the dependence of the exchange-correlation (xc) 
energy E^c on the current density jpo-- This dependence is known to cause only a small effect 
in the electronic ground state of quantum dots. The xc vector potential then drops out, and 
the Kohn-Sham equation as derived in Ref. |TT] becomes 

f ie 

+ n'^^^^f' + ^exto(r) + v^{v) + t;,ea(r) 
= ej^i)j„{v) , (2) 

where as usual Wexto- and Aext denote the external scalar and vector potentials and fH(r) and 
^xco-(r) are the Hartree and xc potentials. In the following, the basic assumption is that the 
quantum dots can be treated as two-dimensional systems. We then use polar coordinates 
and write the external potential as a sum of the bare potential of the quantum dot plus the 
Zeeman term: 

Ve^tair, V?) = VAot{.r, V^) + ^ g*fiBBa . (3) 
Here, fiB = eh/2m, a = ±1, and the specific form of Vdot{r,ip) will be discussed in section 




II Q below. In turn, the external vector potential is given by 
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Aext(r,¥?) = — r , (4) 

corresponding to a uniform magnetic field B perpendicular to the dot. Next, we expand the 
Kohn-Sham orbitals in polar coordinates as 

^A^) = T.Mr)e-'^\ (5) 

where the (pjiair) are real. The spin density is then given by 

n^(r) =n^{r,!f) = J2 fja\^jA^)\'^ 
j 

= E/.-<^.'-(O0.7v(r)e^('-'>, (6) 

jW 

where the fjo- denote occupation numbers of the orbitals. In the following, all calculations 
will be performed at a small but finite temperature, T = 4.2 K, in order to avoid convergence 
problems at small magnetic fields. The occupation numbers are then given by the thermal 
distribution 

and the chemical potential is fixed through the relation J2ja fja = ^- Next, the confining 
potential of the dot is expanded as 

v^ot{r,^) = J2^<^ot,i{r)e-'''' . (8) 

Similar expansions are made for the Hartree and xc potential. For the dot potential, the 
angular components 



1 ^'^ 



dipe''^Vdot{r,'p) (9) 



in general have to be obtained through straightforward numerical integration, and similarly 
for the xc potential (for the latter we use the local-density approximation in the parametriza- 
tion of Tanatar and Ceperley0). For the Hartree potential, one finds 

oo 

vn,iir) = 27re*2 j ^'^r' 0,n<.(r')0,n'.(r') hir, r') , (10) 

C inn' r\ 



n — n'—l 



where e* is the effective electronic charge, and the Ii{r,r') involve integrals over Bessel 
functions: 

oo 

hir, r') = Jdq Ji{qr)Ji{qr') . (11) 
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The Ii{r,r') can be expressed in terms of hypergeometric functions that can be further 
reduced to complete elliptic integrals (see Ref. |^ for the case / = 1). Inserting everything 
into the Kohn-Sham equation, we finally obtain 




im*c^ ' 2 

+ [^dot,fc(r) + Vu,k{r) + V^caA^) 4>jl-ka{r) 
k 

= eja 4>jUr) . (12) 

Eq. (|T2|) couples the angular components of the Kohn-Sham orbitals (^. In practice, we of 
course limit the expansions of the orbitals and potentials to a finite number of components, 
—L, . . . ,1, . . . , L, where L of the order 5 is sufficient to give convergence for the cases under 
study. Solution of the Kohn-Sham equation (|12[) is then accomplished by discretization on 
a logarithmic radial mesh with N^j.^^ ~ 100 grid points. 



B. Linear response 

The FIR absorption spectra as measured in Refs. are proportional to the photoab- 
sorption cross section 

(7{uj) = lma{uj) , (13) 

where the dipole polarizability with respect to, say, the x-axis 

2e f 

a{uj) = — — / (fr X ni{r,u!) (14) 



is obtained from the linear density response ni(r, u) of the quantum dots to an external field 
of the form 

^i(r,c^) = ^— , (15) 

where Eq is the amplitude of the electric field strength. In order to calculate ni(r,co') = 
X^cr '^i(T(r, Col), we have to solve the linear spin-density response equation, which reads as 
follows: 

n.(r,.)^/.Vx.(r,r',.){.,.(r',.) 
+ E / d'r" [y^^ + /xc..(r', r", u;)^ n.^r", u;)^ 

(16) 

i.e. we have to solve two coupled integral equations for ni| and rii^. The Kohn-Sham 
response function 
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Xa{r,r,uj) = 2^{fka - fja) ; — 

7^ Ska- Bja + UJ + tr] 

is diagonal in the spins. Let us now expand the density response as 



N 



imp 



(17) 



(18) 



n=~N 



where in practice of course N < L. Inserting the form (|^) of the Kohn-Sham orbitals, the 
response function is set up in the following way: 



N 



ilip —il'ip' 



Ll'=-N 



where we define 



with 



Xll'a[r,r ,UJ) - 2^[fka - fja)- 

j^k ^ku — Bjcr + UJ + 11] 



^jkiir) = Y <Pji{r)<Pkv{r) 
1,1' 

l'-l=L 



(19) 



(20) 



(21) 



For the xc kernel we use the adiabatic local-density approximation (ALDA) 



/xc,ar(r',r",u;) 



dua drir 



6{r' - r") 



no(r') 



N 



6{r'-r") Y Un.'ar{r')e'"'^ 

m'=-N 



(22) 



where Cxc is the xc energy density of the homogeneous two-dimensional electron gasjli and 



27r 



fxc,m'(TT (^ ) 



dip 



d^e 



271 J dUa- drij 



(23) 



no(r') 



In the ALDA, the xc kernel is frequency independent and real (for a recent discussion of 
alternative expressions for /xc see Ref. ^). The imaginary part of the Kohn-Sham response 
function (^) thus has to be put in by hand. In the following, we choose a value of = 0.1 
meV, corresponding to about 0.1-1% of the excitation energies under study. 
Inserting everything into the response equation (|16]), we obtain 

y N 

nii„{r,uj)= dr'r Y Xii'a{r,r' ,uj) vw^ir' ,uj) 
1'=-^ 

OO OO j\j 

+ An\*^Yhr'r' jdr"r" Y XiVa{r,r' ,uj) m,r{r\uj)k,{r\r") 

r i i V=-N 

°9 N 

+ 27rY J dr'r' Y Xii'<T{r,r' ,uj) riimrir' ,uj)Uc,{i'-m)aTir') . (24) 



l'm=-N 
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Solving for the density response ni(r,co') = ni|(r,a;) + nii{r,u) for a given value of u thus 
requires inversion of a complex matrix of dimension 2A^grid(2A^ + 1) = 1400 (for = 3 and 
-^grid = 100), which poses no problem in practice. 



C. External potential 

Let us now turn to the specific form of the bare confining potential of the quantum dot, 
Vdot{r,ip), used to construct the electronic ground state in Sec. [11 A| . We first of all restrict 
ourselves to considering only potentials that have inversion symmetry. In other words, the 
expansion (H) of Vdotir,ip) contains only angular components Udot,/('") with even /. We can 
then replace e*'"^ in eq. by cos{lip). Furthermore, to reduce the computational effort we 
shall restrict the values of / to / = 0, ±2, ±4. 

We first consider the circularly symmetric part of the confining potential, fdot,o(^)- This 
component describes the degree of non-parabolicity of our quantum dot. Pfannkuche and 
Gerhardsi assume a form 



Vdotfl{r) 



m 



/ 2 2 



(25) 



where Ui -C ujq. A different approach has been chosen by Broido et alB: they construct 
'f^dot.o('") as the electrostatic potential associated with a two-dimensional jellium disk of radius 
R and uniform positive areal charge density Their result is 

^^dot,o('") = ^0 - 4:e*'^n+RE , r < R 
Vdot,o{r) = Vo- Ae*'^n+r 



E 



R 



1 




r > R. 



(26) 



where vq = 27re*^n+ R. Here, K and E denote complete elliptic integrals of the first and 



second kind. It is easy to see that 
coefficients as 



51) is the small-r expansion of 



m*R 
3ne*^n+ 

16m* R^ 



identifying the 



(27) 



The main difference between the two forms of fdot,o(^) is that (p5|) grows as for large r, 
whereas (^) approaches the constant vq as 1/r. This difference is of less importance for 
the electronic ground state, since the two potentials are very similar in the interior region of 
the dot where the electronic density is concentrated, but it can be expected to substantially 
affect the electronic excitations. 

We now turn to the components of fdot that break the circular symmetry. In Refs. 

this is accomplished by including terms of square symmetry, i.e. proportional to x^y^. 
Similarly, one can add on terms proportional to or y"^, describing elliptic elongation of 
the dot along the x or y axis.0 In this manner, one arrives at 
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^'dot,±2(r) = ar^ 

VdotMr) =br\ (28) 

which introduces two more adjustable parameters a and b, in addition to c^o and uJi. 

Again, an alternative approach to the construction of fdot(^;V') is to start out with a 
flat jellium disk of uniform positive charge n^, but this time with a non-circular shape 
parametrized as R{f)- Once a particular form for R{(p) has been chosen, the associated 
electrostatic potential is calculated as 



e*\+ I dip' f , (29) 

1^2 _|_ _ 2rr' cos{ip — (p') 







where now 

= e*'^n+ dip' R{ip') . (30) 
Jo 

The r'-integral in (^) can be performed analytically, and the remaining integration over (p' 
has to be done numerically for general R{p>'). If R{^') = const., one recovers the previous 
result (PBD- The next step is then to construct the angular components fdot,/('^) using eq. 
(§), which in general requires a second straightforward numerical integration. 
The large-r behavior of the dot potential (EUI) is found to be 



Vdotir ^ oo,ip) = VQ ^ / d(p'R{(p') , (31) 







independent of ip. This means that for large distances only the / = component of fdot('", f) 
survives (approaching the constant Vq), and the higher-Z components go to zero. This is again 
in contrast with the form (|28|) for the non-circular components of the dot potential, which 
(unphysically) keep increasing with distance. As noted before in the case of the circular dot, 
this difference is not expected to have a large impact on the electronic ground state, but 
it may become important for higher excitations. We therefore conclude that in general it 
is preferable to work with dot potentials constructed according to eq. (^Of ), thus avoiding 
effects caused by an unphysical behavior in the large-r region. 



III. RESULTS AND DISCUSSION 

The experiments presented in Refs. were performed on self-assembled InAs quantum 
dots embedded into GaAs. The diameter of the dots is estimated to be about 20 nm and the 
height to be about 7 nm. As outlined above, we treat the quantum dots as two-dimensional 
systems. Within our model, we also ignore the presence of the wetting layer. 

From their measurements, the authors of Refs. deduce an effective mass m* = O.OSme 
(where rUe is the bare electronic mass). Furthermore, we take the effective charge as e* = 
e/y/e (e is the bare electronic charge), using e = 15.15 for the dielectric constant, i.e. the 
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bulk value of InAs, and we employ an effective (^-factor g* = —0.44. For the curvature of 
the bare confining potential of the dot close to its center, we take a value of uq = 45 meV, 
which leads to n+ = 0.7 x 10^^ cm"^ via relation (^) (for R = 100 A). 

The specific form of the bare confining potential of the dot can now be constructed using 
one of the two approaches discussed in Sec. |L1 C| . In the following, our choice is to construct 
■ydotl?", V^) as the electrostatic potential associated with a disk whose radius is parametrized 
as 

R{(p) = Rq + R2 sin^if + R4 sin^(y9 cos^yj . (32) 

For i?o, we take the estimated radius of the dots, i.e. Ro = 100 A. The parameter R2 
indicates an elliptic elongation of the quantum dot along the j/-axis. We adopt a value of 
i?2 = 5 A, as estimated In Ref. |[ For the parameter i?4 that causes an anisotropy of square 
symmetry, no direct experimental numbers are available. In the following we choose i?4 = 83 
A, so that the value of R4 sin^y? cos^ip is at most 20% of Rq. The resulting shape of the dot, 
a rectangle with rounded-off corners, is shown as inset in Fig. ^ For comparison, a circle 
with radius Rq is also indicated. 

In Fig. we plot the angular components Vdot,i{r) of the bare confining potential of 
the quantum disk parametrized by (|32|). The top part shows the circularly symmetric part 
Vdot,o{i"), together with a parabohc potential that would correspond to the case fdot('") = 
m*ujQr'^ /2. The bottom part shows the / = 2 and I = 4 components. Note that if the 
elongation of the dot is along the x-axis [replace i?2sin^v' by R2Cos'^ip in eq. (|32D], then 
Vdot,2('") changes sign. From Fig. |1] it is evident that the deviation from circular symmetry 
affects the confining potential mainly in the region around the edge of the dot, whereas the 
inner region of the dot sees a nearly parabolic potential. 

Let us now discuss our main numerical results. Fig. Q shows the calculated peak positions 
of the photoabsorption cross section cr(to') versus applied magnetic field, for a quantum dot 
with N = 2 electrons. Here and in the following, we assume the quantum dots to be 
elongated along the x-axis. We then calculate the photoabsorption spectra for two different 
polarizations of the FIR radiation, in x- and y-direction, respectively. The symbols in Fig. 
^ denote the average of the two spectra. For comparison, the full lines show the expected 
behavior of lj± according to eq. (|I|), with a fitted value of ujq = 46.8 meV. We see that for 
higher magnetic fields, the calculated peaks follow the simple law of eq. (|^). For small B, 
however, some deviations occur, and a splitting of about 2 meV remains even for i? = T, 
in accordance with observation.ii 

For N = 2, the quantum dot contains a full s-shell,0 and the behavior is very similar to a 
parabolic dot. As soon as the p-shell gets occupied, however, deviations from the parabolic 
case become much more pronounced. In Fig. ^ we plot the FIR peak positions for an = 3 
quantum dot versus B. The intensities of the absorption peaks are approximately indicated 
by the size of the symbols. The two top figures show results for the non-circular dot, for 
two different polarizations of the FIR radiation: parallel to the direction of elongation of the 
dot (the X-axis) and perpendicular to it. In both cases we find very rich spectra. For small 
magnetic fields, there are substantial differences between the two polarizations, for B > 6 
T these differences disappear. 

The bottom part of Fig. |^ shows the spectra for a circular (but still non-parabolic) dot. 
These results have been obtained by setting Vdot,i{r) = for / = ±2, ±4, but using the same 
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'i^dot,o('^) as in the two figures above. By comparison between the lower and the two upper 
parts of Fig. we can now clearly distinguish between those effects related to breaking the 
circular symmetry and those caused by non-parabolicity. The latter leads to a splitting of the 
ui^ mode into two almost equally strong subpeaks, 0;+^ and a;+\ separated by about 5 meV, 
plus the appearance of a weaker signal connected with the uj- mode and approximately 10 
meV below it. We mention that similar results have been previously obtained by Hawrylak 
and coworkers. 00 

(3) 

By introducing a non-circular anisotropy, the a;+ mode acquires a third subpeak uj\ in 
agreement with the experimentally observed behavior (see Fig. 3 of Ref. |^). We find that the 
evolution of this signal with magnetic field is different from that of the two other subpeaks 
of ct;+: it becomes weaker with increasing B and its separation from and is growing. 

To understand the origin of the three- fold splitting of u!+, it is helpful to resort to a 
simple single-particle picture. Fig. ^ shows the energy levels of a single electron confined in 
the bare potential of our quantum dot at = 12 T, versus angular momentum quantum 
number / (see footnote 0). The dashed lines connect energy levels to which one can assign 
the same principal quantum number n. The distribution of energy levels shown in Fig. ^ is 
very similar to that for a parabolic dot, which is governed by the formula 




E^i = {2n+\l\ + l)Juji + ^ - ^ (33) 



(see the discussion in the review article by AshooriEl). In contrast with the strictly parabolic 
case, the vertical distances between the levels in Fig. ^ are not constant. 

In this simple picture (which remains qualitatively valid for > 1 even if the electron- 
electron interaction is included), the subpeaks of can be identified with single-particle 
transitions. For circular symmetry, these are governed by the selection rule for dipole transi- 
tions, Al = ±1. In this simplified scenario, 0;+'' and arise from Is lp_ and lp+ 2s 
transitions, as indicated in Fig. 11. Breaking the circular symmetry means that the selection 
rule can be violated: we find that the uj\_ mode originates from the transition lp_|_ ld_, 
i.e. Al = —3. The oscillator strength of this mode is of course directly related to the degree 
of anisotropy, which in our case is only small. 

We have found that the position of u)^^ with respect to and uj^^ is insensitive to 



small variations in the choice of the parameters R2 and R4 in eq. (32). A weak signal at 
uj\^ will be present even if only one of the two is nonzero. However, the value of i?2 = 5 
A is more or less dictated by the experimentally observed splitting between 0;+ and U- for 
A^ = 2 at zero magnetic field, see above. One then observes that choosing a finite value for 
i?4 increases the intensity of the third subpeak. 

The central result of this work is presented in Fig. ^. It shows a comparison between the 
calculated and measured! peak positions in the photoabsorption spectra for quantum dots 
with 1 < A^ < 6 electrons at i? = 12 T. We find that the main experimental features are 
reproduced by the calculation. For A^ = 1,2, the system behaves very similarly to a circular 
parabolic dot, as noted before, i.e. there are only two signals at uj^ and uj^. As soon as 
the p-shell becomes occupied, i.e. for A^ > 3, the u+ mode splits up. As explained above, 
the uj^^ signal (indicated here by the open circles) is related to a breaking of the circular 
symmetry. 
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We note that the calculation yields a splitting between the three subpeaks of uj+ that is 
greater than the one found in the experiment. Also, for = 6 the experiment yields only a 
single signal at ujj^. These differences between theory and experiment are to be attributed 
to the simplified nature of our model that treats the self-assembled dots as two-dimensional 
disks. In particular, the observed differences at = 6 are most likely due to our neglecting 
the presence of the wetting layer. 

The main effect of the wetting layer is to introduce a continuum of states above a certain 
energy thresholdi, limiting the number of bound states localized in the self- assembled dot. 
As the number of electrons grows, the energy levels are shifted towards higher energies due to 
the increasing interaction energy, and more and more states are pushed into the continuum, 
up to a point where no additional electron can be bound. From the absence of uj\ and 
uj\ for A^ = 6 observed in the experiment, we infer that the relevant states involved in 
the transitions (2s and ld_) would in reality fall into the wetting layer continuum and the 
transitions would lose most of their strength. However, for A^ < 6 these states must still 
be located in the discrete part of the energy spectrum, since the associated transitions are 
present in the experiment. 

This effect has been accounted for in Ref. by using a truncated basis of only few bound 
states in the numerical diagonalization, and it was indeed found that for A^ = 6 there is only 
a single peak at a;+\ However, the calculations in Ref. |1^ assumed circular symmetry and 
produced at most a twofold sphtting of 0;+ for A^ = 3, 4, 5, in contrast with experiment. We 
therefore conclude that, in spite of the deviations from experiment mentioned above, our 
calculations clearly establish that the presence of u;^'' for A^ = 3, 4, 5 is due to a breaking of 
circular symmetry of the quantum dots. 



IV. CONCLUSION 

In this work, we have developed a theoretical description of collective charge-density 
excitations of non-circular quantum dots in a magnetic field. The computational scheme 
presented here allows one to obtain information on the geometry of quantum dots from 
their electronic response properties. 

Our specific aim was to reproduce and explain recent experimental results, obtained 
for self-assembled InAs quantum dots. In these FIR photoabsorption measurements, one 
detects a threefold splitting of the upper branch ojj^ of the collective charge-density mode 
in a magnetic field, depending on the number of electrons present in the dot. Our study 
has shown that these experimental features are closely related to the shape of the dots: 
in addition to non-parabolicity (due to the finite dot radius), it is essential to account for 
anisotropy effects, leading to a non-circular confining potential. To our knowledge, we have 
presented here the first fully self-consistent spin- density-functional calculations of ground 
state plus linear response for anisotropic quantum dots with up to 6 electrons. 

With our calculations we were able to reproduce and explain the main features of the 
FIR spectroscopy measurements conducted in Refs. ||,^. However, the agreement was not 
fully quantitative. This may be attributed to the simplicity of our model, which treats 
self-assembled, lens-shaped (with a possible elliptic or pyramidal distortion) quantum dots 
sitting on a wetting layer as two-dimensional quantum disks. It may safely be expected 
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that this model leads to much better quantitative results for quantum dot devices that are 
produced by mesa-etching techniques, also known as vertical quantum dotJll. 
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FIGURES 




FIG. 1. Angular components Vdot,iif) of the bare confining potential of a quantum disk 
parametrized by eq. (|32|). The shape of the disk is shown in the inset, together with a circle of 
radius Ro = 100 A. Top: / = component (full line). For comparison, the dashed line shows a a 
parabolic potential with the same curvature in the interior region. Bottom: / = 2 (full line) and 
/ = 4 (dashed line) components. 
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FIG. 2. Calculated peak positions of the photoabsorption spectrum for a non-circular, 
non-parabolic quantum dot with N = 2 electrons. The full lines show the results one would 
obtain for a circular parabolic dot, see eq. (||). 



15 



1 r 



1 r 



1 r 



CO 



70 - 



% 60 



% 60 



> 60 



50 



40 



30 



X 



50 - 



40 - 



30 - 



(2) 

+ 



H h 



70 - 



Y 



50 - 



40 - 



30 



• • • • • • 



H h 



70 - 



C 



: i ; t : : 



_j I I I I I I 

2 4 6 8 10 12 

B{T) 



FIG. 3. Calculated peak positions of the photoabsorption spectrum for a quantum dot with 

N = 3 electrons. Bottom: non-parabolic but circular dot. Top and middle: non-parabolic, 
non-circular dot, with polarization of the FIR radiation parallel (X) and perpendicular (Y) to 
the direction of elongation. The mode splits up into three subpeaks, lo^\ lo^^ and lo^\ as 
indicated. 
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FIG. 4. Energy levels of a single electron in the bare potential of the quantum dot from Fig. 
1^ at = 12 T. The levels are drawn versus their angular momentum quantum number I, and 
dashed lines connect levels with the same principal quantum number n. The arrows indicate the 
single-particle transitions that constitute the three- fold splitting of uj^. 
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FIG. 5. Comparison between the calculated and experimental peak positions in the photoab- 
sorption spectra for quantum dots with 1 < < 6 electrons at = 12 T. The open circles indicate 
those signals that we find to be related to a broken circular symmetry of the dots. 
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